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Spin Chirality in a Molecular Dysprosium Triangle: the Archetype of the 

Non-Collinear Ising Model 
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Single crystal magnetic studies combined with a theoretical analysis show that cancellation of 
the magnetic moments in the trinuclear Dy'^"'" cluster [Dy3(^3-OH)2L3Cl(H20)5]Cl3, resulting in 
a non-magnetic ground doublet, originates from the non-coUinearity of the single ion easy axes of 
magnetization of the Dy^"*" ions that lie in the plane of the triangle at 120° one from each other. 
This gives rise to a peculiar chiral nature of the ground non-magnetic doublet and to slow relaxation 
of the magnetization with abrupt accelerations at the crossings of the discrete energy levels. 



PACS numbers: 71.79.Ej, 75.10.Jm,75.50.Xx, 75.30.Cr 



Molecular nanomagnetism has provided benchmark 
systems to investigate new and fascinating phenomena 
in magnetism[lj, y] like magnetic memory at the molecu- 
lar level[3|, quantum tunneling of the niagnetization[J, 5|, 
or destructive interferences in the tunneling pathways 6| . 
In this field rare earth ions like dysprosium(III) are 
currently investigated because of their large magnetic 
anisotropy and high magnetic moment [7|. In the course 
of our synthetic efforts to obtain new molecular nanomag- 
nets based on rare-earth ions we recently obtained the 
trinuclear Dy3+ cluster [Dy3(/X3-OH)2LaCl(H20)5]Cl3 
(where L is the anion of ortho- vanillin) Q , hence ab- 
breviated as Dya, which possesses an almost trigonal 
(C3/1) symmetry (see Figure [T] and EPAPS for more 
information) Q . 

Preliminary powder magnetic measurements measure- 





FIG. 1: (color online) Left:View of the molecular structure 
of Dy3 cluster where the hydrogen atoms, the chloride coun- 
teranions and the solvent molecules of crystallization have 
been omitted. Right: Schematic view of the spin structure of 
the Dy3 triangular cluster and of the local easy axes orienta- 
tion in respect of the laboratory XYZ reference frame. 



ments on two compounds containing such Dy3 clusters, 
but differing in their crystal packing and intermolecular 
contacts, revealed an identical magnetic behavior where 
the magnetization vs. field curve at low temperature 
is almost flat, suggesting a non magnetic ground state, 
but suddenly increases to its saturation value, Ms, at H 
= 8 kOeQ- This is in contrast to what commonly ob- 
served for antiferromagnetic triangular clusters that show 
a multi-step magnetization curve [lO[. A possible expla- 
nation is that the Dy'^"'' ions are characterized by almost 
ideal Ising anisotropy with the single ion easy axes lying 
in the plane of the triangle at 120° from one another, as 
represented in the scheme of Figure [T] To the best of our 
knowledge an experimental realization of this simple but 
fascinating spin structure is unprecedented. 

To verify our hypothesis larger crystals (of size ca. 1 
mm'^) of one of the two compounds were grown according 
to @] through very slow evaporation of the solvent. This 
allowed an accurate face indexing of the crystal on the 
X-ray diffractometer and the investigation of the mag- 
netic anisotropy by using an horizontal sample rotator 
in the SQUID magnetometer (see EPAPS for experimen- 
tal details) [01 . Scans in different crystallographic planes 
allowed us to determine the three magnetic anisotropy 
axes, denoted as X, Y and Z. The two structurally 
equivalent Dys molecules in the unit cell have the Dy3 
planes almost perpendicular to Z and one side of the tri- 
angle parallel to Y (see Fig. [T]). Magnetization vs. field 
curves along these axes are given in Figured^. Along X 
and Y a sudden jump around 8 kOe is observed while an 
almost linear but weaker magnetization is observed along 
Z. In Figure[2b the temperature dependence of the mag- 
netization along the three axes measured at 1 kOe, and 
thus before the jump to saturation, are shown. The in- 
plane X and Y directions are very similar and M tends 
to zero at low temperatures, confirming a non-magnetic 
ground state, while a weaker signal is observed along Z. 
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The observed behaviour has been modelled using the 
canonical formalism of the statistical thermodynamics 
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FIG. 2: (color online) a) Field dependence of the magnetiza- 
tion measured along the X, Y, and Z axes at r=1.9 K. In 
the inset a magnified view of the high field region is reported, 
b) Temperature dependence for the susceptibility along X, Y, 
and Z axes with an applied field of 1 kOe. The solid lines 
represent the calculated value with Eq. (2) and the best fit 
parameters discussed in the text. 



FIG. 3: (color online) a) Zeeman splitting of the levels due 
to the application of the field along the X axis (bisector) 
calculated with the Spin Hamiltonian (1) and the parameters 
indicated in the text, assuming = 90°. b) Same diagram 
calculated with (j) = 0° . The spin structure for each state is 
schematized by the arrows. The non-degenerate states are 
highlighted by dash lines. 



and taking into account that T^y^^ ions have a very large 
magnetic single-ion anisotropy due to the crystal field 
splitting of the ^-ffi5/2 ground state. In a first approx- 
imation each Dy^+ ion, which is supposed to have the 
doublet ground state well separated in energy from the 
other excited Stark sublevels, can be represented by an 
effective spin 5 = 1/2 and the Dys system can be mod- 
elled by an Ising Hamiltonian: 
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where for each Dy^+ ion the Zi local axis is consid- 
ered to be in the plane of the triangle at 120° one from 
each other as schematized in Figure [1] The basis set of 
the full system consists therefore of 8 vectors: | tTT)i 

I itt), I Tit), I TTi), I Mi) , I iti), I lit), I iii), 

where, however, up and down refers to the local z axes. 
Actually, these vectors are already the eigenvectors of 
(1). The best-fit parameters of the power data were 
jzz/ks = 10.6(4)7^ and g^ = 20.7(1). This last value 
confirms that the ground doublet for the Dy'^+ ions is well 
described by \J = 15/2, mj = ±15/2). In fact J = 15/2 
results from the coupling of L = 5 and S = 5/2 and 
thus gj = 4/3, which, for mj = ±15/2, gives an effective 
gyromagnetic factor of 20. 



It is interesting to focus on the energy of the eight 
eigenstates of the Hamiltonian (1) as a function of the 
applied magnetic field in the Dya plane, as shown in Fig- 
ure [3l At zero field two states with a zero net magnetic 
moment, are degenerate with energy equal to — 3J2z/4, 
whereas the other 6 states are also degenerate with en- 
ergy equal to jzz/4. In an applied magnetic field the 
energy of the last six levels depends on the angle the 
magnetic field forms with the local easy z axes. Let us 
suppose that the field is applied along one bisector of 
the triangle as shown in Figure [3l Two limiting scenar- 
ios can be observed depending if the local easy axes are 
along the edges of the triangle (Fig. [3^) or along the bi- 
sectors (Fig. EJd). a jump to saturation magnetization at 
the first level crossing. Hi, is expected in both cases, as 
indeed observed in Dya. In general, thanks to the struc- 
tural non-collinearity of the easy axes, the two states of 
the ground doublet have opposite vortex chirality, with 
clockwise or anti-clockwise rotation of the spins. 

This simple model cannot however reproduce the weak 
linear increase in the magnetization above the step, 
visible in the inset of Figure [5^. We have therefore 
formulated a new Hamiltonian model by considering 
for each Dy'^"*" that the ground state doublet is |J = 
15/2, mj = ±15/2) and that the first excited doublet 
is|J = 15/2, mj = ±13/2). These states can be effi- 
ciently admixed by the J^ and J_ operators of the Zee- 



man interaction. To avoid ovcr-paramcterization we have 
limited the treatment to the first excited doublet. The 
expression for the new Hamiltonian is: 
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where a,; runs over the local axes of the i Dy ion (a; ^ 
Xi, Hi, Zi). The first term of (2) comes from considering 
an isotropic (Heisenberg) exchange between two Dy ions, 
where {di,dk) is the angle between the ai and ak local 
axes. The last term describes the single-ion anisotropy 
and d is the zero field splitting between \J = 15/2, mj — 
±15/2) and I J = 15/2, mj = ±13/2) states of each Dy3+ 
ion. The best fit, again performed on the powder data, 
provided.? = -0.092(2)ii', g = 1.35(1) and 6 = 102(5)i5:. 
The g value is now close to 4/3 as expected for J — 
15/2. The large separation between the two Kramers 
doublets is in good agreement with what is reported in 
the literature for Dy'^+ [ll|, [l2| . To reproduce the single 
crystal data the angle 4> between the X axis and the 
local z anisotropy was allowed to vary freely. The best 
simulation is obtained with the angles 4> = ±17(1)° ± 
ri60° (n=0,l,2 ..), where the periodicity results from the 
symmetry of the cluster. We have also evaluated the 
intra-molecular dipolar contribution to j as a function of 
(j), and in the most favourable configuration {(/) = 90°) it 
can only account for about half of the observed value. 

Susceptibility measurements using standard induction 
coils in alternating magnetic field allowed us to estimate 
the relaxation rate from the frequency dependence of the 
imaginary component of the susceptibility x" assuming 
that, according to the Debye model, at the maximum 
of x" '^s CO curve the simple relation r = l/w is valid 
(see EPAPS for more details) [1|. When the ac field is 
applied in the plane of the triangle the relaxation rate 
decreases on lowering the temperature following an Ar- 
rhenius law, r = tq exp(A/fcsr), with to=2.5(5)x10^^ 
and A =36(2) K. At temperatures below ca. 7 K the 
relaxation increases less rapidly, as shown in Figure [31 

The Arrhcnius behavior of the relaxation time char- 
acterizes a class of molecular materials called Single 
Molecule Magnets, SMMs[i3|, where the easy axis mag- 
netic anisotropy generates a barrier for the reversal of 
the magnetization giving rise at low temperature to mag- 
netic bistability and memory effect of pure molecular 
origin [l|, 13, Sj. Dys, however, represents the first example 
where such a slowing down of the magnetization dynam- 
ics occurs even if the overall magnetization lies in an al- 
most isotropic plane rather than along an easy axis. That 
is because the system can be better schematized by three 
interacting SMMs. In the past pairs of weakly coupled 
Mn4 clusters have been widely investigat ed Il4l ilSj and 
demonstrated to show quantum coherence |16|. However 
a fascinating new situation is observed in Dys because 
the two states of the ground doublet are distinguished 
by a different spin chirality and they cannot be related 



one to each other by a simple exchange of the magnetic 
sites. The last is verified even in the special non-chiral 
case of = 0°. 

One of the striking aspect of SMMs is that relax- 
ation at low temperature can occur through an under- 
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FIG. 4: (color online) Arrhenius plot of the temperature de- 
pendence of the relaxation time obtained from the frequency 
dependence of the ac magnetic susceptibility in zero and ap- 
plied static field far from a level crossing. The lines represent 
the linear fit according to the Arrhenius law of the highest 
temperature data. 




5 10 

H (kOe) 



FIG. 5: (color online) a) Field dependence of the real com- 
ponent at 5K of the ac susceptibility measured with the field 
along X at logarithmic spaced frequencies. The empty circles 
represent the experimental static susceptibility obtained by 
derivation of the experimental magnetization curve, b) Field 
dependence of the relaxation time measured at T=7 K with 
the field applied along Y . In the inset the Zeeman splitting 
calculated for H parallel to Y with the best fit parameters 
discussed in the text. 



barrier mechanism jl[. The tunneling is particularly effi- 
cient close to zero field, where the maximum degeneracy 
of the states is observed [J, [5[. The apphcation of a static 
field has indeed a strong influence on the dynamics of 
the magnetization of Dys, as shown in Fig. [5^ where 
the real component of the susceptibility, %', for different 
frequencies in 1-1000 Hz range is shown as a function of 
the field. Around zero field and the first level crossing. 
Hi, the dynamic susceptibility approaches the equilib- 
rium value, i.e. that obtained from the derivative of the 
static magnetization curve recorded at the same temper- 
ature, while for intermediate fields a strong frequency de- 
pendent reduction, arising from the impossibility to fol- 
low the oscillating field, is observed. Similar results have 
been obtained along Y, as shown in Fig [5b where the 
field dependence of the relaxation rate at r=7 K for this 
orientation is given. The acceleration of the relaxation 
occurs at a field slightly smaller than that corresponding 
to the maximum of the static susceptibility, suggesting 
that only the first crossing is relevant for the dynamics 
(inset of Figure [Sja) . Far from the level crossings, i. e. at 
if =3 kOe, the relaxation times, given in Fig 51 show a 
significant increase of the barrier to A=120(4) K. 

The present results are well rationalized by the model 
we have previously developed. The response in the ac 
field mainly involves transfer of the population from/to 
the ground doublet to/from the first magnetic excited 
state, implying the reversal of one spin inside the triangle. 
Far from any level crossings this seems to occur through 
an Orbach process involving the first excited Kramers 
doublet of each of the Dy'^+ ions, and its energy gap 6, 
estimated from the static properties, is found indeed to 
be in good agreement with the observed activation energy 
A. On the contrary, at H=Q the reduction of the barrier 



suggests that an alternative mechanism takes place. Tun- 
nelling of the magnetization has already been observed 
for other lanthanide-based SMM and is attributed to the 
admixture in zero external field of the ground Kramers 
doublet states, made possible by the hyperfine interac- 
tion [ITI, lis] ■ At the first level crossing the situation of 
zero local field is reestablished for a Dy'^+ ion resulting 
in the observed fast dynamics. 

To conclude, Dya has revealed to be a benchmark sys- 
tem to investigate non-coUinearity in Ising systems. It 
is worth stressing that non-coUinearity is a key feature 
of molecular magnetism, where complex building-blocks 
characterized by low symmetry are assembled in more 
symmetric architectures. Dy3 combines the slow dynam- 
ics of Single Molecule Magnets and the level crossing ob- 
served in antiferromagnetic rings [19l |2(3J . Both types of 
systems are currently being investigated for their poten- 
tial application in quantum computation[2l|, [22, [23, [241 
and systems with a non-magnetic nature of the ground 
doublet state could be used in order to reduce decoher- 
ence effects due to the fluctuation of local magnetic fields. 
In the ideal case when the spins lie exactly on the plane 
of the triangle the dynamics involving the ground dou- 
blet is however not directly accessible with magnctom- 
etry and requires to be further investigated with more 
sophisticated techniques, for instance using local probes 
like muons or neutrons. 
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